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Abstract 

Energy transfer from photoexcited excitons localized in silicon nanoparticles to adsorbed oxygen molecules excites 
them to the reactive singlet spin state. This process has been studied experimentally as a function of nanoparticle size 
and applied external magnetic field as a test of the accepted understanding of this process in terms of the exchange 
coupling between the nano-Si exciton and the adsorbed O2 molecules. 
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Background 

Since the discovery that photoexcited siUcon nanoparti- 
cles can act as energy donors to molecular oxygen accep- 
tors and can thereby excite oxygen to a highly reactive 
singlet state [1-3], there has been much work on the 
potential exploitation of this process. Applications that 
have been demonstrated range from photodynamic can- 
cer therapy [4,5] to optically activated reactors in chemical 
engineering [6]. 

In early work, it was demonstrated that the efficiency 
of the energy transfer process is sensitive to an exter- 
nally applied magnetic field [2] (the energy transfer effi- 
ciency may be monitored by its quenching of the nano-Si 
residual photoluminescence), and this provided key evi- 
dence for the understanding of the process as a result of 
exchange coupling between an exciton confined within 
a silicon nanoparticle and an adsorbed oxygen molecule 
(the Dexter exchange mechanism). The applied magnetic 
field B lifts the spin degeneracy of both the exciton and 
oxygen spin manifolds; both oxygen molecules and silicon 
excitons will then relax predominantly into their lowest 
energy spin states at temperatures T for which gfiB^ > 
kT where g = 2.0 is the gyromagnetic ratio and /x^ is 
the Bohr magneton. The energy transfer process between 
these lowest energy spin states has a low probability due 
to angular momentum selection rules, so that the effect 
of the magnetic field at low temperatures is to suppress 
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the energy transfer from the exciton to the molecular oxy- 
gen. As a result, the silicon photoluminescence intensity 
is restored towards the intensity observed when oxygen is 
not present. 

Although earlier investigations proposed this model [2], 
the response to a magnetic field has not been investigated 
or modelled quantitatively in terms of the dynamics of 
the energy transfer and other excitation and relaxation 
processes. Furthermore, the dependence of the efficiency 
of the process on oxygen concentration has never been 
investigated. Here, we show results of experimental inves- 
tigations at lower oxygen concentrations than used previ- 
ously, and we set out a preliminary model which makes 
some simplifying assumptions but which has the features 
required to describe our experimental data. This model 
is a starting point for a full theoretical description of 
the energy transfer phenomenon and can be expanded 
to model the energy transfer process as a function of, 
for example, nanoparticle size. Even at the present level 
of approximation, the modelling turns out to be a fairly 
complicated task requiring a large set of input param- 
eters, though many of these are available in the litera- 
ture; some we use have been estimated as part of the 
present work. 

Methods 

The samples were produced in the form of porous silicon 
layers (thickness of approximately 8 \im) on bulk crys- 
talline substrates by conventional electrochemical etching 
from wafers consisting typically of p-type boron-doped 
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CZ < 100 > silicon with resistivities of 1 to 25 ^2 cm. 
Room temperature anodization was performed in a 1:1 
solution of 49% aqueous HF and hydrous ethanol; the 
porosity p was varied by variation of the current (10 
to 40 mA/cm^) and was determined by fitting of the 
Fabry-Perot interference fringes in a broad-band optical 
reflectance measurement [7] to be typically p = 63% 
to 70%. The etched layers were left attached to the sub- 
strates for better mechanical strength and were glued to a 
copper cold finger with heater and thermometer resistors 
attached. The samples were held either in a continuous- 
flow cryostat (base temperature of approximately 10 K) 
or a superconducting magnet in superfluid helium (base 
temperature of approximately 1.5 K). The magnetic field 
was varied up to 6 T and was oriented either parallel 
or perpendicular to the sample normal. The orientation 
of the field plays no role in the following experiments, 
in which the optical polarisation of the photolumines- 
cence (PL) emission was not analysed. The effects we 
discuss here depend only on the magnitude of the induced 
Zeeman splittings in the exciton and oxygen triplet states 
(polarisation-dependent studies are under way at present). 
In both cryostats, the cold finger could be raised to the top 
of the cryostat to expose the cold sample briefly to oxy- 
gen gas and it could be heated whilst in vacuum to desorb 
oyxgen. PL was excited by a continuous wave solid state 
diode laser (wavelength approximately 450 nm, power 
approximately 5 mW at the sample, with a weakly focused 
laser spot, size a few hundred microns) and detected with 
an intensified CCD camera and compact single-grating 
spectrometer. 

Results and discussion 

Four typical PL spectra at 1.5 K for a porous silicon sam- 
ple exposed to a low oxygen concentration are shown in 
Figure 1 (spectra were recorded at 0.5-T intervals, but for 
clarity, we omit the spectra at intermediate fields). The 
broad luminescence band corresponding to a wide dis- 
tribution of silicon nanoparticle (NP) sizes is observed 
[8-10]; this band is similar in shape to that obtained in 
the absence of oxygen but is lower in intensity. The over- 
all intensity of the PL band increases by about 20% as 
the applied magnetic field is increased to around 4 T 
and then ceases to increase further. This behaviour differs 
quite markedly from the first reported experiments using 
a magnetic field, where the oxygen concentration was high 
enough that PL above the threshold energy of 1.63 eV 
for singlet oxygen production was still completely sup- 
pressed even at fields as high as 10 T and the field-induced 
recovery of the PL intensity was only observed below 
1.63 eV [2]. 

Figure 2 shows the PL spectra obtained at higher oxygen 
concentrations (Figure 2) in a second piece of the porous 
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Figure 1 Photoluminescence of porous silicon containing a low 
concentration of molecular oxygen. Photoluminescence (PL) 
spectra of a porous silicon sample exposed to a small quantity of 
oxygen gas are shown for magnetic fields of 0 to 6 T. The sample was 
held in superfluid helium at 1 .5 K, and the PL was excited with 
450-nm (2.76 eV) continuous wave excitation. The vertical dashed line 
indicates the threshold energy, above which photoexcited excitons in 
the silicon nanoparticles have sufficient energy to excite the 
adsorbed oxygen from its triplet to its singlet ^ E state. 



silicon sample used to obtain the results of Figure 1. It 
is not possible to measure quantitatively the oxygen con- 
centration adsorbed on the silicon NPs, but the much 
stronger quenching of the PL gives a clear indication that 
the concentration is higher than in the case of Figure 1. 

There are two notable features: Firstly, the strongest 
quenching of the PL occurs precisely for NPs having an 
exciton energy equal to the oxygen to transition 
energy of 1.63 eV. Secondly, the spectra show a large 
number of other sharp downward-pointing peaks or dips 
which originate from the enhanced energy transfer to oxy- 
gen for NPs whose exciton energies differ from 1.63 eV by 
energies corresponding to one or more momentum- and 
energy-conserving phonons (located at K and F points 
of the silicon phonon dispersion, respectively). These 
phonon effects have been discussed elsewhere, where 
details of the relevant phonon energies are given [3]. Two 
prominent dips of this type can be seen near 1.9 and 
2.0 eV; these are also related to energy transfer to oxy- 
gen but will be discussed in future work; here, we shall 
model only the energy transfer process without phonon 
participation. 
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Figure 2 Photoluminescence of porous silicon containing a high 
concentration of molecular oxygen. Photoluminescence (PL) 
spectra of a porous silicon sample exposed to a larger quantity of 
oxygen gas than in Figure 1 are shown for magnetic fields of 0 to 6T. 
As in Figure 1 , the sample was held in superfluid helium at 1 .5 K, and 
the PL was excited with 450-nm (2.76 eV) continuous wave excitation. 
The vertical dashed line again indicates the threshold energy for 
energy transfer, at which the quenching of the PL is particularly 
efficient. Other structures arise from energy transfer processes in 
which phonons participate. 



Figure 2 demonstrates that significant PL is again 
observed above the threshold for energy transfer to oxy- 
gen, even at this higher oxygen concentration. Further- 
more, the PL both above and below this threshold shows a 
much stronger recovery of intensity as the magnetic field 
is increased, by factor of about 3 times, and unlike the case 
of Figure 1, the recovery of the PL has not saturated up to 
a magnetic field of 6 T. 

The differences between Figures 1 and 2 point to an 
interplay between the rates for the physical processes 
(light absorption, radiative recombination, spin relax- 
ation, and energy transfer) that control the shape of the 
PL spectrum. These processes are indicated schematically 
in Figure 3, which serves as a guide to the rate equation 
model we develop below. Figure 3 summarises the situa- 
tion of NPs with oxygen present, for which there are four 
possible states (represented by the four boxes): the oxygen 
molecule can be in either a singlet or a triplet state, and the 
NP may or may not contain an exciton. Optical pumping 
creates excitons, whilst PL emission and energy transfer 



processes annihilate them. Only energy transfer generates 
singlet oxygen, whilst spin relaxation (or infrared PL) pro- 
cesses return the oxygen to the triplet ground state. In 
the rate equation model for these processes, the photoex- 
cited populations of the separate spin states of the excitons 
and the oxygen molecules are treated explicitly, taking into 
account the spin dependence of the energy transfer to O2, 
the radiative exciton recombination rate, the processes of 
thermal excitation and spin-lattice relaxation that lead to 
population redistribution between the spin states for a 
given silicon NP, and the rates of relaxation from singlet to 
triplet oxygen states. 

Silicon nanoparticles without oxygen 

At the low measurement temperatures necessary for 
magneto-optical experiments (we use 1.5 K), we know that 
oxygen is not able to desorb from the nanoparticle sur- 
faces (experimentally, we have to heat the sample to about 
80 K before oxygen is released and can be pumped away, 
at which point the PL intensity recovers completely). We 
can therefore divide the NPs into two separate popula- 
tions: those which are in contact with oxygen (represented 
in Figure 3) and those which are not. We write the propor- 
tion of NPs which do not have adsorbed oxygen molecules 
and which do not currently contain an exciton as hq; exci- 
tons are created in these in one of the three triplet exciton 
states (index / = 1 ... 3) with equal pumping rates P/3 to 
generate fractional populations Ui, The photoexcited NPs 
can de-populate only by radiative emission with rates ro, ri 
for mj = 0, my = ±1, respectively (note that, here, we set 
these equal; we will consider the consequences of these 
being different in a future work), spin-lattice relaxation to 
spin states lower in energy (y^y), or thermal excitation to 
spin states higher in energy by A// (/// = yexp(—Aij/kT)), 
Note that A/y is dependent on the magnetic field since 
it arises from the Zeeman splitting of the exciton states; 
this leads to a magnetic field dependence of yij. Non- 
radiative relaxation processes may also contribute to the 
triplet exciton relaxation at low temperatures [11] but 
would enter into our model in the same way as the radia- 
tive decay rates and so are not included explicitly. Under 
these assumptions, the steady state solution of the rate 
equations for the fractional populations Ui, hq yields the 
following result (Equation 1): 

iP/3)no - uiin + 712 + 713) + U2y2i + ^^3/31 = 0 
(P/3)no + uiYu - U2(ro + 721 + K23) + ^^3K32 = 0 
(P/3)no + uiyi3 + i/2K23 - min + K31 + 732) = 0 
no -\- ui -\- U2 -\- U3 = I - F, 

(1) 

where F is the total fraction of NPs with adsorbed oxygen. 
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Figure 3 Schematic overview of energy transfer from photoexcited excitons in silicon nanoparticles to absorbed oxygen molecules. 

Optical excitation (green arrows, 'pump') generates excitons confined in silicon nanoparticles that can recombine to emit photoluminescence (red 
arrows, 'PL') or can transfer energy to those absorbed oxygen molecules that are in the triplet ground state (black arrow, 'energy transfer'). Excited 
oxygen molecules in the singlet state can return to their ground state (blue arrows, 'relaxation') via emission of luminescence and/or non-radiative 
relaxation processes. 



Silicon nanoparticles with oxygen 

We now consider the second population of NPs, those 
which are in contact with oxygen. We write the propor- 
tions of NPs which do not contain an exciton as nj, where 
; runs over the three possible oxygen triplet states. As 
above, excitons are created in these NPs in one of the 
three triplet exciton states (index / = 1 ... 3) with equal 
pumping rates P/3 to generate fractional coupled exciton- 
oxygen populations riij. The exciton radiative recombina- 
tion and spin-lattice relaxation terms are as above, and we 
introduce a spin-lattice relaxation and thermal excitation 
term between the oxygen triplet states analogous to yij 
(Pij), Note, again, that pij is in general a function of mag- 
netic field and depends on both zero-field and Zeeman 
terms (shown in Figure 4). We must also account for NPs 
in which the oxygen is in the singlet state and no exciton 




Figure 4 Energy level diagram for the energy transfer from 
photoexcited silicon nanoparticles to oxygen molecules. Left: the 
triplet (bottom) and singlet (top) levels of molecular oxygen in a 
magnetic field, showing the zero-field splitting between the mj = 0 
and xUemj = ±1 levels; right: the ground state (bottom) and triplet 
exciton (top) states of a silicon nanoparticle in a magnetic field. In this 
model, the exchange coupling between the two triplet states results 
in energy transfer between three possible pairs of states, of which 
one example is indicated by the vertical and curved arrows. 



is present (the condition of an NP after energy transfer 
and before relaxation of the oxygen, with population He) 
and NPs in which an exciton has been excited whilst the 
oxygen is still in the singlet state (populations Wj), 

Finally, we introduce the energy transfer process which 
is the focus of this work through the rate tij. In the simplest 
approximation, as represented in Figure 4, the magnetic 
field and the principal axis of the oxygen molecule can 
be taken to be parallel; to model the behaviour with a 
random distribution of angles between these directions is 
substantially more complicated (requiring an average over 
the relative orientations and a calculation of the mixing of 
spin states) and will be discussed in future work. Here, our 
aim is to investigate what can be achieved with a realis- 
tic set of parameters in a comparatively simple model. The 
matrix tij here has the following form in order to impose 
the overall conservation of spin angular momentum, 
Am/ = 0: 



(2) 



As in the previous subsection, we present the steady 
state solutions of the resulting 15 rate equations plus the 
condition that the total number of NPs with adsorbed 
oxygen remains constant. The first sets of expressions 
(Equations 3 to 5) represent the generation and loss of 
excitons in NPs with adsorbed triplet oxygen; the exis- 
tence of two triplet entities gives nine possible joint spin 
states, so that nine equations are required. 

(P/3)ni - (ri + tn + yu + yis + ^12 + Pis) ^ii 

+ y21^21 + YSimi + ^^21^12 + fti^is + (R/3)wi = 0 

(P/3)ni - (ro + t2i + 721 + ¥23 + P12 + Pis) ^21 

+ Yunil + y32W31 + ^^21^22 + Al^23 + {R/3)w2 = 0 

{P/3)ni - (ri + + 731 + y?>2 + P12 + ^^13) ^31 

+ Yunii + 723^21 + hin?>2 + Psims + (^/3)m/3 = 0 




(3) 
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(P/3)n2 - (n + ti2 + yi2 + yi3 + P21 + «i2 

+ 721^22 + 731^32 + Punn + ^2^13 + W3)wi = 0 
(P/3)f72 - (ro + ^22 + K21 + K23 + ^^21 + ^^23) ^22 

+ yunu + y32W32 + ^^12^21 + ^^32^23 + (^/3)W2 = 0 

(P/3)n2 - (n + ^32 + K31 + K32 + ki + P23) n32 

+ n3Wl2 + K23W22 + ^^12^31 + feW33 + (^/3)W3 = 0 

(4) 

(P/3)W3 - (ri + ti3 + yi2 + yi3 + fti + ^2) «13 

+ 721^23 + 731^33 + Pisnn + yS23Wl2 + (R/3)Wi = 0 

(P/3)W3 - (ro + % + K21 + K23 + Psi + P32) n23 

+ m^l3 + y32W33 + ^^13^21 + ^3^22 + (^/3)W2 = 0 

(P/3)W3 - (ri + ^33 + Ysi + 732 + ^31 + ^^32) ms 

+ n3'^13 + m'^23 + )^13'^31 + yS23W32 + W3)W3 = 0 

(5) 

The next set of equations (Equation 6) represents the 
optical pumping and de-excitation of NPs with adsorbed 
oxygen in its singlet state; the three equations arise from 
the three exciton states. 

iP/3)ne - (ri + + yi3 + 712) wi + y2iW2 + ysiws = 0 
(P/3)ne + 7121^^1 - (ro + + 723 + Yii) ^2 + ¥32^3 = 0 
(P/3)ne + XlsWi + ¥23^2 - (n+R+ Y31 + 732) VI/3 = 0 

(6) 

The final set of equations represents the generation and 
loss of NPs with triplet oxygen but no exciton; the rate 
R expresses the oxygen relaxation from singlet to triplet 
state. 

{R/3)ne -{P + P12 + Pu) ni + ^21^2 + y^3iW3 

+ nnii + ri^si + roW2i = 0 
{R/3)ne + Puni - (P + ^21 + fe) «2 + ^2^3 

+ nnu + ri«32 + ron22 = 0 
iR/3)ne + ^13^1 + few2 -(P + P31 + P32) n3 

+ ri«i3 + ri«33 + ro«23 = 0 

As stated above, the remaining equation (Equation 8) 
imposes the requirement that the total fraction of NPs 
with adsorbed oxygen should remain constant at F, With 
this condition, we have a fully determined system and can 
solve for all 16 variables in this equation. 

«e + «1 + «2 + «3 + 1^1 + 1^2 + 1^3 

+ mi + ni2 + ni3 + W21 + n22 (8) 

+ ^23 + «31 + n32 + «33 = F 



We can sum all the exciton radiative processes in order 
to obtain an expression for the PL intensity /pl as follows: 

/pL = n (^13 + ^33) + ron23 

+ n (ni2 + + ^0^22 

+ n {nil + n3i) + ro«2i (9) 
+ n {wi + W3) + rQW2 
+ n (ui + U3) + roM2 

and this expression can be evaluated as a function of mag- 
netic field; note that riij, Wi and, in principle, Ui are all 
functions of magnetic field through the field dependence 
of Yij and fiij. 

Comparison to experiment 

The above model does not account for phonon-assisted 
processes and therefore is strictly only valid for NPs emit- 
ting PL at the threshold energy of 1.63 eV. In fact, this 
is not a serious limitation, since the degree of recovery 
of the PL in a magnetic field is similar over a PL energy 
range wide in comparison to a phonon energy. It is beyond 
the scope of this work to discuss the energy dependence 
of the transfer process in detail, and so we extract only 
the PL intensities at 1.63 eV from the spectra of Figures 1 
and 2 and plot them in Eigure 5 as a function of magnetic 
field, normalized to the PL intensity at zero field. This 
normalization eliminates the difficulties associated with 
considering absolute PL intensities and will facilitate the 
comparison of data from different samples. 

Eigure 5 also shows calculated results based on the 
above model, in which we take a set of parameters based 
on the recent literature. These are summarised in Table 1. 
Eor the two sets of experimental data, we maintain all 
parameters at the same values, except for those associated 
with the energy transfer process itself: these are F, which 
expresses the proportion of NPs without oxygen, and the 
transfer rate t, which decreases as the probability of an NP 
having multiple O2 molecules available increases. 

The fraction F of NPs with adsorbed oxygen was varied 
from 0.75 (Eigures 1 and 5, blue) to 0.85 (Eigures 2 and 5, 
red), and 1/t varied from 10~^ to 10~^ s. More work is 
needed before we would attempt to interpret these param- 
eters directly, but we note that these transfer times are 
in good agreement with previously measured values [12], 
and as is necessary for the evenly matched competition 
between radiative recombination and energy transfer, they 
are comparable to the radiative lifetimes l/ri,l/ro [13]. 
In the simulations, we also varied the temperature, since 
the field at which the PL recovery approaches saturation 
is sensitive to the relationship between gfiB^ and kT. As 
can be seen from Eigure 5, the simulations agree well with 
the experimental results taking the nominal experimental 
temperature of 1.5 K. We will report elsewhere on studies 
of the excitation intensity dependence of the effect; there, 
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Figure 5 Comparison of experimental data and results of the 
rate equation model. Solid points: tine ratio of tine PL intensity at 
magnetic field 1(B) to that at zero field l(B = 0) (red circles and blue 
squares: high and low O2 concentrations, respectively); lines: 
predictions of the rate equation model for l(B)/l(B = 0) keeping all 
parameters constant except those related to the oxygen 
concentration and for a series of temperatures (upper to lower 
curves) of 1 .5 to 4.5 Kin 1 -K steps. 



Table 1 Parameters used in modelling (inverse rates, in 
seconds) 



This work Typical Source 





Low O2 


High O2 






Silicon NP 












10-5 


10-5 


10-5 to 10-2 


[13] 




10-5 


10-5 








10-7 


10-7 








1/45 


1/45 






Oxygen 










F 


0.75 


0.85 








4 X 10-3 
2 X 10-7 


4x10-3 
2x10-7 








10-5 


2x10-7 


2.6 X 10-^ 


[12] 



we find we must take into account an increase in temper- 
ature for high excitation intensities (here, these were the 
same for Figures 1 and 2 and were low). 

Conclusions 

Using the simple model set out above, the dependence 
of the photoluminescence spectra of silicon nanoparti- 
cles with adsorbed oxygen molecules has been studied 
and it is shown that a realistic set of parameters can 
give an adequate description of the recovery of the PL 
intensity with increasing magnetic field, confirming the 
proposed spin-dependent exchange-coupled mechanism 
for the energy transfer process. In particular, one set of 
parameters can describe the behaviour of the magnetic 
field dependence for high and low oxygen coverage of the 
sample by changing only the parameters directly relevant 
to the energy transfer process. This represents the first 
detailed and quantitative investigation of magnetic field 
effects in the photogeneration of singlet oxygen by use of 
silicon nanoparticles and provides a model which can eas- 
ily be expanded in order to investigate the dependence of 
the energy transfer process on nanoparticle size, excita- 
tion intensity, and temperature; this work is in progress. 
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